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ABSTRACT 

Aims. We present a cosmic shear analysis and data validation of 15 square degree high-quality 7?-band data of the Garching-Bonn 
Deep Survey obtained with the Wide Field Imager of the MPG/ESO 2.2m telescope. 

Methods. We measure the two-point shear correlation functions to calculate the aperture mass dispersion. Both statistics are 
used to perform the data quality control. Combining the cosmic shear signal with a photometric redshift distribution of a galaxy 
sub-sample obtained from two square degree of UBVRI-hand observations of the Deep Public Survey we determine constraints 
for the matter density flm, the mass power spectrum normalisation erg and the dark energy density 57a in the magnitude 
interval R £ [21.5,24.5]. In this magnitude interval the effective number density of source galaxies is n = 12.5 arcmin"^, and 
their mean redshift is z = 0.78. To estimate the posterior likelihood we employ the Monte Carlo Markov Chain method. 
Results. Using the aperture mass dispersion we obtain for the mass power spectrum normalisation a"g = 0.80 ± 0.10 (la 
statistical error) at a fixed matter density f2m = 0.30 assuming a fiat universe with negligible baryon content and marginalising 
over the Hubble parameter and the uncertainties in the fitted redshift distribution. 
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1. Introduction 

Measuring the weak gravitational lensing effect induced 
by the tidal gravitational field of the large-scale struc- 
ture (LSS) of the Universe, also called cosmic shear, is 
a powerful way to explore statistical properties of the 
LSS and, thus, cosmological models. In contrast to other 
methods cosmic shear directly probes the total matter 
distribution of the Universe independent of the baryonic 
content. Due to the need for high quality data, it has 
only been in recent years that groups have begun to use 
this method. The first cosmic shear measurements were 
published in the beginnin g of 2000 llBacon et alJl2000t 
Van Waerbeke et al.ll200ollKaiser et alJl2000HMaoli et alJ 
200lt IWittman et al.ll2000|) . Although some constraints 
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on cosmological parameters were obtained from these 
measurements, their statistical errors, due to the small 
areas of the surveys (^Ideg^), and systematic errors 
were large. More recently, published cosmic shear mea- 
surements use larger survey areas (~ lOdeg'^) and/or 
deeper fields to obtain smaller statistical errors and ad- 



the cosmic shear signal fe.g. 


Van Waerbeke et al 2601I 


l2002ll200f4lR,hodes et alJ200ll 


2004HRefreEier et al.l2002l: 


'Hamana et al.' '2003'; 'Bacon et al.' '2003*; 'Hoekstra etal 


2002; 


Brown et al. 2003; Masscv et aX, .2005; .Jarvis et al. 


^oosl 


l200ilHevmans et al.ll200f5ll . 



We are now entering a new phase of wide-field sur- 
veys which are wider and/or deeper so that systematic er- 
rors may begin to dominate over statistical errors. Many 
groups have been working on the problem of reducing sys- 
tematic errors in cosmic shear measurements and have 
presen ted its impact o n cos mological parameter e stima- 
tions. iHoekstral l)2004l) and IVan Waerbeke et alJ l)2005l) 
discussed the influence of an imperfect PSF-anisotropy 
correction on the cosmic shear signal. The impact of 
the redshift distribution sampling error has recently 
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been s tudied bv lMa et all l)2006|) and lVan Waerbeke et alJ 
l)2006|) . A further, more fundamental, source of bias is the 
correlation between the weak gravitational shear of dis- 
tant galaxies and the intrinsic shape of foreground galaxies 
(e.g. iHirata fc Seliak l'2004':'H evmans et al.l2006 bY In ad- 
dition, the so-called Shear TEsting Programme (STEP), a 
world-wide collaborative project to improve the accuracy 
and reliability of all weak lensing measure ments has been 
initiat ed in mid 2004 (see first results in iHevmans et alJ 

l2nnfia|). 

In the present paper we first present the results from 
a simple test of our pipeline using synthetic images. 
Specifically we give details about the changes we have 
made to it due to the lessons we have learnt so far from 
STEP and present the influence of the changes on the 
cosmic shear signal. To test the pipeline we used 3 deg^ of 
synthetic images obtaine d by ray tracing through A^-body 
simulations (.Hetterscheidt et al A cosmic shear 

analysis is performed using the E- and B-mode decom- 
position of the aperture mass (Map-) statistics. We show 
that within the uncertainties we recover the cosmic shear 
signal and that our pipeline does not create artificial B- 
modes. 

We then report on the results of a cosmic shear analysis 
of the Garching-Bonn Deep Survey (hereafter: GaBoDS) 
data set, obtained with the Wide Field Imager (WFI) of 
the MPG/ESO 2.2m telescope on La Silla, Chile. The 
GaBoDS data set consists of 18.6 deg^ high quality R- 
band observations. Most of the data were obtained in the 
framework of a virtual survey. That is, we compiled the 
data set from the ESO archive. The rest of the data set 
was obtained by our GO programmes. 

We effectively use 15deg^ of the total data set to cre- 
ate catalogues of galaxies for the shear measurements. 
Furthermore, the data set includes 2 dcg^ of deep UBVRI- 
band observations from the Deep Public Survey (hence- 
forth: DPS). This yields photometric redshift information 
for 8% of all GaBoDS lensing objects in the magnitude 
interval R G [21.5,24.5]. The DPS lensing galaxies with 
photometric redshift information are used to estimate the 
redshift distribution which in turn is considered for the 
cosmic shear analysis. 

To test the reliability of our measurements we perform 
four basic tests on the systematics in our data. In addi- 
tion, we calculate the ellipticity cross-correlation between 
uncorrected stars and corrected galaxies to quantify the 
residuals of the PSF anisotropy correction. Furthermore, 
we utilise the E- and B- mode decomposition of the Map- 
statistics to check for possible systematics in the data. 

The signal of the aperture mass statistics, the shear 
two-point correlation functions and their covariance ma- 
trix are calculated for the cosmic shear analysis. In con- 
trast to other cosmic shear analyses, we estimate the mass 
power spectrum normalisation, erg, using the Monte Carlo 
Markov Chain (MCMC) technique. 

The paper is organised as follows. In Sect.|21we briefly 
describe the theory of cosmic shear. We also discuss the 
statistics that we use to constrain cosmological parameters 



as well as a control of possible systematics in the data. The 
data and the weak lensing catalogue creation are presented 
in Sect. 01 The analysis of synthetic images is presented 
in Sect. ^ In Sect. we elucidate in detail the different 
tests on the systematics in our data. In Sect. El we give 
the results from our cosmic shear analysis and present 
the cosmological parameters calculated using the MCMC 
technique. In addition we discuss several possible sources 
of systematic errors. In Sect. [3 we give a summary and 
outlook. 

2. Cosmic Shear 

In the following we describe the relevant quantities for 
cosmic shear analyses. We use t he standard lensing nota- 
tion t hrou ghout the paper See iBartelmann fc SchneideJ 
l|200l|) and lSchneider et alJ l)2006|) for a more detailed in- 
troduction. 

2.1. Relation between power spectrum and shear 
estimators 

The power spectrum P„ of the projected density field (con- 
vergence k) is defined as 



^.(0 = 



4c4 



a2(w)^ \fK{w) 



I 



(1) 



where Hq is the Hubble constant, c the speed of light, 
w is the comoving radial coordinate, a is the cosmic scale 
factor, P3D is the three-dimensional mass power spectrum, 
I is the modulus of the two-dimensional wave vector, and g 
is the weighting function accounting for the source redshift 
distribution of galaxies Pzdz ~ p^dw, 



9{w) 



™« „ Uw'-w) 
/k(w') 



(2) 



see for instance iKaised l)l998l) or ISchneider et al.l l)l998j) . 

The quantity fxiw) is the comoving angular diameter dis- 
tance and wn is the distance to the horizon. 

In this work we use the shear two-point correlation 
functions, and the aperture mass dispersion, 

(Mfp), to estimate cosmological parameters. For com- 
pleteness we also give the shear dispersion, (7^). All quan- 
tities are linearly related to the convergence power spec- 
trum by 



27r Jo 



{ML){0o) = — dllP,{l)W,piWo) 



2tt 



and 



{i^){0o) = 7r diip,{i)WTM), 

2,71 



(3) 
(4) 
(5) 

(6) 
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with 

Wi^{v)^Jo{v),w^Av)^Mv), (7) 

Wth(?7)=4^ and W,p{rj) ^ 576^^, (8) 

where J„ denote the n"^-order Bessel function of the first 
kind. 

Th e aperture mass, Map, was introduced bv lSchneideil 
and is defined as the spatially filtered projected 
mass distribution, k, inside a circular aperture of angular 
radius at a position ^, 



A4p(0 = J d'9K{e)u{\e-c\), 



(9) 



where J7 is a radially symmetric continuous compensated 
weight function given iui^chncidcr ct al. (1998). The aper- 
ture mass Map can also be expressed in terms of the ob- 
servable tangential shear jt 



Map(C) = I d'9jt{9;C)Q{\e-C\), 



(10) 



where Q is a poly nomial filter function given in 
ISchneider eraD l)l998|) . 

The cosmic shear signal can be decomposed into 
a curl-free (E-mode) a nd a curl (B-mode) component. 
ICrittenden et all l)2002f) showed that B- modes do not con- 
tribute to Map and E-modes do not contribute to AI± 
which is defined as 



Mi(C) - / d207x(0;OQ(|e-CI), 



(11) 



where 7x is the cross component of the shear. As the aper- 
ture mass unambiguously separates E- and B-modes, it 
can therefore be used to reveal possible systematics in the 
data. 

The aperture mass dispersion can directly be obtained 
from the observed images by placing, for example, a reg- 
ular grid on the data field. For each grid point the Map- 
and Afj^-values can be measured and then their dispersion 
calculated. However, with this method it is impossible to 
obtain unbiased estimates for (Map) and {Mj_) because of 
gaps, holes and borders in the data fields. We therefore 
calculate the dispersions directly in terms of the shear 
two-point correlation functions, which are directly related 
to the observable quantities 7t and 7x : 



?±(^) = (7t7t)(^)±(7x7x)(0), 
?x(e) = (7t7x)(e). 



(12) 
(13) 



Due to parity invariance the correlation function is 
expected to vanish. 

The aperture mass dispersion is related to the shear 
two-point correlation functions, by 



(M^)(^o)4/f 



(14) 



(Mi)(.o)4/f 



Another quantity defined by 



(15) 



(16) 



can also be used as an indicator for systematics as it is 
unaffected by cosmi c shear. The functions T± can be ex- 
pressed analytically ( Schneider et alJl2'002l) and vanish by 



definition for 9 > 29^. Thus (Mfp), (M^) and {Ml) can 
be obtained directly from the correlation functions ^± in 
a finite interval. In Eq. 1)16(1 we choose r_ because it has 
no change of sign in contrast to T+ . 

In this work we utilise the aperture m ass dispersion to 
separ ate E- and B-modes (as suggested in ISchneider et alJ 
Throughout this paper we use (Mfp) and E inter- 
changeably as well as {M\) and B. 



2.2. Shear estimators in practice 

The estimators of the two-point shear-shear correlation 
functions in practice are 



X;f,- WiWj (et,iet J ± ex ,iex J ) {9) 



i±{0) - 

with 
A,,(0) 



1, e - M/2<\9,~9j\<e M/2 
0, otherwise. 



(17) 



where et and ex are the tangential and cross components 
of the corrected galaxy ellipticities and Wi is the weighting 
factor (details of the correction and weighting scheme can 
be found in Sect. 13.31 and i3.4|l . All pairs of background 
galaxies within a distance oi 9 — A0/2 and 9 + A6'/2 are 
considered, where A0 is the angular width of the bin. The 
practical estimators of the aperture mass dispersion are 



M 



e+(0„)T+P^ +|_(0„)T 



Mx 



„=i 

2 m 



(18) 



with the centre of the bins 9n — {n — 1/2) A9 and the 
aperture radius 9o. 

3. Data and galaxy catalogues 

3.1. The GaBoDS data set 

The Garching-Bonn Deep Survey (GaBoDS) data set 
was obtained with the Wide Field Imager (WFI) of the 
MPG/ESO 2.2m telescope on La Silla, Chile. The cam- 
era consists of eight 2k x 4k CCDs with a pixel size of 
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15 /im, corresponding to a pixel scale of 0'.'238 in the sky. 
The field-of-view is 34' x 33'. However, due to the applied 
dither pattern the effective field-of-view can be as large as 
40' X 40'. 

The total data set consists of 18.6 deg^ of high qual- 
ity i?-band observations. About 80% of the data were ob- 
tained by a virtual survey (ASTROVIRTEL project^-^), 
meaning that the data set was collected from the ESO 
archive. For data mining the archive, the search engine 
querator'^ was developed by Pierfederici ( 2001) to address 
the needs of astronomers looking for multicolour imaging 
data. Fields were selected by four criteria: (1) the expo- 
sure time was at least 3.5 ks, (2) the seeing was better than 
1", (3) the field was not affected by bright foreground ob- 
jects and (4) neither the field nor adjacent fields contained 
previously known very massive structures. For further de- 
tails of the data min ing of the ESO archive we refer to 
ISchirmer et alJ ((20031) . The rest of the data set was ob- 
tained by our GO programmes (tableau. 

The data set is heterogeneous with respect to the 
Vega-limiting magnitude in the i?-band (ranging between 
25.0 mag and 26.5 mag; 5 a sky level measured in a circular 
aperture of 2" radius), seeing (ranging between 0'.'75 and 
1'.'20), observing strategy (dithering/no dithering) and 
co-addition process (different programmes were used for 
co-addition: ElSdrizzle and swarp'', see table The 
positions of the fields are randomly distributed on the 
southern hemisphere at high galactic latitudes, however 
some of the fields form small patches in the sky. The 
survey is separated into six data sets mainly depending 
on the data source, see table and lA.ll Note that the 
ESQ Distant Clus ter Survey CEnisCS. IClowe et ablEooi 
IWhite et al.l2005l) data set contains distant galaxy cluster 
candidates (z = 0.5. ..0.8) that have been identified using 
the red cluster sequence. Since we perform our lensing 
analysis with those objects lying in the magnitude inter- 
val R G [21.5,24.5], we are presumably not sensitive to 
those high redshift clusters. We assume in the following 
that these fields do not bias our cosmic shear survey sig- 
nificantly towards high-density regions. 

Each of the WFI images observed with dither pattern 
(w 70% of all fields) cover an area of about 0.35 deg^. 
However, for the cosmic shear analysis the edges are 
trimmed off (with a low effect on the total number den- 
sity of background sources) resulting in an effective area 
of 0.25deg2 per field. 

Furthermore, the data set includes 2deg^ of 
deep VB V7?/-bai id ob servations from the DPS 
(|Hildebrandt et alJ l2006|) . which yields photometric 



^ ASTROVIRTEL cycle 2: Erben et al., Gravitational lensing 
studies in randomly distributed, high galactic latitude fields 
^ http: / / www.stecf.org/ astrovirtel 
^ http:/ /archive. eso. org/querator 

* lErben et al.l lj200^ showed that there is no significant dif- 
ference in the size of the PSF or the PSF anisotropy pattern 
using ElSdrizzle or swarp. 



Table 1. Data code, depth of the fields (time t of all co- 
added exposures), total field area, effective number density 
noff and data source (see table lA.II for details of individual 
fields). Key: *: co-addition with swarp, ^: single exposures 
are not dithered, "'': fields are rejected for the final cosmic 
shear analysis. 



code 


depth t 


area 


Woff 


Data source 




[ks] 


[deg2] 


['-'] 




COMBO 


10<t<56 


1.25 


16.5 


COMBO-17 


OWN 


7.0<t<22 


3.25 


12.0 


own obs., 
MPE IR group, 
ESO archive 


DEEP* 


3.9<t<9.3 


3.25 


11.5 


EIS + 
own obs. 


Bg!" 


4.2<t<7.5 


2.5 


11.5 


ASTROVIRTEL 


C0+'* 


4.0<t<4.8 


2.0 


10.0 


ASTROVIRTEL 


CL 


t = 3.6 


3.0 


12.0 


EDisCS 



redshift information for 8% of the objects considered for 
the cosmic shear analysis. 

3.2. Data reduction 

The data reduction has been performed with a nearly fully 
automated, stand-alone pipeline which has been developed 
by our group to reduce optical and near-infrared images, 
especially those taken with multi-chip camera s. A detailed 
description of the pipeline can be found in lErben et alJ 
f2005). As weak gravitational lensing was the main science 
driver, the pipeline algorithms are optimised to produce 
deep co-added mosaics from individual exposures obtained 
from empty field observations. Special care has been taken 
to achieve an accurate astrometry to reduce possible ar- 
tificial PSF patterns in the final co-added images. Since 
keeping the PSF anisotropics small is mandatory for cos- 
mic shear measurements, individual exposures are rejected 
from the co-addition process if one or more CCDs exhibit 
an anisotropy larger than 6 % in either of the ellipticity 
components. More information a bout the exact d ata re- 
duction process can be found in ISchirme^ ()2004|) . Note 
that for the image co-addition of the DPS data set (code 
name in our work: DEEP, table the programme swarp 
was used in contrast to the rest of the GaBoDS data set 
that was co-added with ElSdrizzle. 

3.3. PSF correction 

The shape of galaxies is influenced by the anisotropic 
PSF. In order to obtain a correct estimate of the she ar 7 
from t he observed ellipticity of galaxies x"*^**,, Kaiser et alJ 
(|l995() developed the so-called KSB algorithm. The algo- 
rithm relates the observed ellipticities x°^^ to the sheared 
source-ellipticities and provides an unbiased estimate of 
the shear. The correction is calculated on the Gaussian 
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kernel weighted second brightness moments Qij 
galaxy with sm'face brightness I{0), 



9,)ie,~e,)iie)w{\e 



-0\ 



with 

w{\e- 



0\ 



2na- 



rexp 



2(72 



of a 



(19) 



(20) 



where the scale length, a, is defined as the half-light ra- 
dius, Th, of the object considered, and the centre is 



9 



_ Jd'^ew{9)ei{e) 
" / (few{9)i{e) ■ 

The ellipticity is defined as 
Qii - Q22 + 2iQi2 



X = 



Qii + Q22 



(21) 



(22) 



Assuming that the intrinsic orientation of galaxies is 
random, the relation between 7 and the weak lens- 

ing regime reads 



(23) 



where is the pre-seeing shear polarisability which de- 
pends on the smear and shear polarisability tensors P*^™ 
and P^^ of the galaxy, and the stellar smear and shear 
polarisability tensors P^"^* and P^^* in the following way 



pg _ pun _ psm^psli*-j — ipi 



(24) 



The tensor PS relates the measured anisotropy-corrected 
galaxy ellipticity to its true value. All quantities are cal- 
culated by means of the observable surface brightness, 
I {6). The quantity q* is the anisotropic part of the PSF 
and is calculated from the raw stellar ellipticity X*i 1* = 

3.4. Catalogue creation 

Raw catalogue. \ye util ise the programme SExtractor 
( Bertin fc Arnout3 Il99fit) to create a raw catalogue of 
objects (source galaxy candidates). The final co-added 
image is first smoothed with a Gaussian kernel of 
2.5 pixel FWHM (FILTER_NAME: gauss_2.5_5x5.conv). 
Then objects are extracted which consist of at 
least N — 5 contiguous pixels (SExtractor parame- 
ter 'DETECT_MINAREA') with a flux greater than 
k — 1.5(7 above the sky level noise (SExtractor parame- 
ter 'DETECT_THRESH'). With these rather conservative 
parameters we reduce the number of noise detections and 
increase the reliability of shape measurements. For the 
object detection the co-added WEIGHT image is taken 
into account as an additional SExtractor argument. It 
provides a full characterisation of the relative noise prop- 
erties for each science image pixel and therefore lowers the 
probability to detect pure noise objects, see Erben ct al. 
i)2005(l . In addition, we set BACK_TYPE=MANUAL and 



BACK_VALUE=0 (our co-added images are sky sub- 
tracted). In this way SExtractor does not model halos 
around bright objects as sky background. It turns out that 
this effectively reduces spurious detections of objects in 
the proximity of bright stars or within their halos. 

Calculation of quadrupole moments, Qij and centroids, 

. The raw catalogue is transfered to the programme 
analyseldac, an adjusted version of N. Kaisers imcat, 
which measures the quadrupole moments of each object 
detected with SExtractor. Due to pixelisation, the contin- 
uous integrals of Eqs. I|19|) and (|21|l are now transformed 
into discrete sums, 

= £ (Q^_Q.)i^Q^_Q.) 1(0) w (|0 - (25) 

and 



Ae^w{\e~e\)i{e) 



(26) 



where 9 is measured, in pixel units, from the source cen- 
troid 6 and Ad — 0.25 pixel (each pixel is divided into 
four sub-pixels). The object centre 6 is determined by it- 
eratively solving Eq. H26I) . where the starting point is the 
SExtractor centroid. The scale in the Gaussian weight- 
ing function, W, is the half-light radius, rh, which is 
equal to the SExtractor parameter FLUX_RADIUS . The 
quadrupole moments, Qij, are calculated from all pix- 
els for which |0 — 0| < 6'niax = 3rh. For non integer pixel 
values of the position 6, the surface brightness I{6i,6j), 
known at pixel positions, is estimated from a linear inter- 
polation over the four nearest pixels to {6i,6j). 

PSF-anisotropy correction. In order to apply Eq. H23|l 
the anisotropy term, q* , has to be known. For the PSF- 
anisotropy correction we therefore utilise stars which are 
point-like and unaffected by lensing. All stars have the 
same half-light radius and therefore show up as a verti- 
cal branch in a magnitude against half-light radius plot. 
Stars which have a brightness of P = 0.5 magnitudes lower 
than saturated stars and which are well above the crowded 
faint magnitude region which contains a mixture of faint 
stars, galaxies and noise detections are selected. Using 
this sample of bright, unsaturated stars, we measure q* 
with a Gaussian filter scale matched to the size o f the 
galaxy image to be corrected. iHoekstra et all lll998|) pro- 
posed this matched PSF-anisotropy correction scheme for 
space-based data. It turned out in STEP 2 (Massey et 
al., in preparation) that this also improves the results for 
ground-based data. 

In the case of the WFI@2.2m instrument the 
PSF anisotropy of the co-added images is rather 
small and varies smoothly over the total field-of-view. 
Discontinuities in the PSF anisotro py are largely absent 
across chip borders (demonstrated in lSchirmer et alJ200,'4 
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ISchirmedEool lErben et alJl2005l) . Therefore we perform 
a second- or third-order two-dimensional polynomial fit 
to q* with 3.5 (T-clipping as a function of position over the 
entire field-of-view. With this fit it is possible to estimate 
the anisotropy kernel, g*, at the position of the galax- 
ies. A basic analysis of PSF anisotropy properties of the 
GaBoDS data set is given in Schirmer et al. (in prepara- 
tion). Note that this anisotropy correction method is not 
perfect as the PSF is only measured at the position of 
stars and a polynomial interpolation is used to estimate 
the PSF anisotropy at the position of galaxies. However, 
our fields comprises several hundreds to some thousands 
of stars used for the anisotropy correction, so that the PSF 
is sampled very well. 

For future cosmic shear analyses aiming at high pre- 
cision constraints for cosmological parameters the PSF 
anisotropy correction method must be improved. Recently 
Jarvis fc Jain (,,2004) have introduced a promising correc- 
tion algorithm, the so-called principal component analysis. 
With this method the dominant PSF anisotropy patterns 
are identified from a large number of exposures and are 
used to model the PSF. In this way, the PSF anisotropy 
correction does no longer depend on the number density 
of stars per field, but by the stacked number density of 
stars across all fields of the survey. 

Calculation of and 7. The tensor is an extremely 
noisy quantity. The diagonal elements of all required ten- 
sors to calculate arc dominant by at least a factor of 
10 compared to the off-diagonal elements, and they are 
approximately equal, so that we can estimate P^ by Pg^l, 
with 

pg = - |tr(P"'>) - tr(P^'»)l , (27) 

2 \ ^ ' tr(P™*) J 

see lErben et all l)200l|) . The stellar smear and shear polar- 
isability tensors P**™* and p^^* calculated for differ- 
ent smoothing scales rj^ and a fifth-order polynomial fit as 
a function of position is performed to tr(P'^™*)/tr(P''^*) 
since it varies over the total field-of-view. As P^ depends 
on these stellar quantities we calculate P^ depending on 
the galaxy siz e, Vh- For P^ we us e the raw, unsmoothed 
values because fErben et al.l l)200l|) pointed out that fitting 
or calculating means does not improve the shear estimates. 

Our final shear estimate is given by 7corr=7//cah 
where 7 is calculated according to Eq. (^5)1 and /cai — 0.88 
is a calibration factor obtained from skymaker and 
shapel et simulations (see STEP 1, 2: iHevmans et al.l 
Massey et al. in preparation). 

Catalogue filtering. All objects for which problems con- 
cerning the determination of shape or centroid position 
occur are rejected (e.g., objects near the border, with neg- 
ative total fiux, with negative Qn + Q22i with negative 
semi-major and/or semi- minor axis, or when the iterative 
centroid determination is not stable after a certain num- 
ber of steps). In addition, we only use those objects with 




^18 20 22 24 26 



R— band [mag] 

Fig. 1. Galaxy counts in i?-band magnitude of the 
GaBoDS for three different exposure time intervals. Solid 
line: deep (10 ks < t < 56 ks), dashed line: medium 
(7ks < t < 10 ks), dotted line: shallow (3.6 ks < t < 7ks). 
All distribution functions are normalised by area between 
18 < P < 24.5. Galaxies with R > 24.5 (shaded area) are 
not considered in this work. 

a half-light radius which is larger than that measured for 
stars and a modulus of the ellipticity (after PSF correc- 
tion) of less than 1.0. 

To obtain a more homogeneous data set, with respect 
to the completeness of the object catalogue of individual 
fields, the raw background galaxy catalogues are compiled 
from objects in the SExtractor isophotal magnitude in- 
terval R € [21.5,24.5]. In this magnitude range we are 
confident to obtain reasonable photometric redshift esti- 
mates for most of the lensing galaxies in the individual 
DPS fields so that we have a reliable redshift distribution 
for all lensed galaxies. In Fig.Qlthe galaxy number counts 
are displayed. The distributions are approximately equal 
below R < 24.5, independent of the exposure time so we 
perform a cut at a magnitude where completeness stops 
for all fields. Hence we would not improve our final cos- 
mic shear analysis significantly by going deeper given that 
faint galaxies are strongly downweighted anyway (Fig. 12). 



Weighting. For each galaxy the next 12 neighbours are 
identified in the magnitude-half-light radius plane and 
the variance, (Jg , of the ellipticity distribution of the sub- 
san iple is calculated (f or more details about this technique 
see lErben et ahlEoOll) . The variance, a^, gives an indica- 
tion of the noise level of this galaxy. From tjg we then 
determine the weighting factor, Wi, as, 



where = 0.16 is the variance of the unlensed galax- 
ies (estimated from the total PSF corrected galaxy sam- 
ple). In Fig. 121 the mean weighting factor is displayed as 
a function of magnitude. In the magnitude interval con- 
sidered, we find that the average relative weight measured 
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PSF correction, same cuts, weighting), see Sect. 13.41 The 
obtained galaxy catalogues are called output catalogues. 

To exclude the effect of false detections we only take 
into account those objects which are present in both, input 
and output catalogues. To obtain a similar number den- 
sity in comparison to the real data we perform a similar 
magnitude cut: mag e [21.5,24.5]. The mean ellipticity 
dispersion of the galaxies of the input and output cata- 
logues is = 0.32 and the mean galaxy number density 
is n = 15 arcmin^^. For a more detailed description of the 
image creati on, galaxy morphology an d m agnitude dis- 
tribut ion see iHetterscheidt et all ((20051) and lErben et alJ 
l|200lh . 



Fig. 2. Average weighting factor calculated from all 
galaxies as a function of magnitude. The weighting fac- 
tors are calculated with Eq. H28() such that galaxies with 
a magnitude of R = 21.5 have a weight of 1. 

bv lHoekstra et ah! l)2005(l is comparable to ours, although 
they use a completely different weighting scheme. 

4. Cosmic shear analysis of synthetic images 

In this section we present the results from qualitative tests 
of our pipeline. For simplicity we divide our pipeline into 
two parts: (1) the part concerning the PSF correction 
and shear estimation and (2) the part which measures the 
two-point correlation function and calculates the aperture 
mass dispersion. In addition, we investigate the impact of 
the changes we have made to part (1) of our pipeline (see 
table 121) on the cosmic shear signal due to the lessons we 
have learnt so far from STEP. 

We performed a cosmic shear analysis of 
twelve synthetic WFI sized images created from 
ray tracing through ACDM iV-body simulations 
{flA = 0.7, no = 0.3, as = 0.9, h = 0.7) using the E- 
and B-mode decomposition of the Map-statistics. These 
simulations were ki ndly made av ailable by Takashi 
Hamana [details see iHamana et al ] §004)]. We created 
twelve initial catalogues of randomly distributed galaxies 
using the programme stuff"" (E. Berlin). The galaxies 
are assumed to be at a fixed redshift z — 1, and are 
sheared according to the shear map of the ray-tracing 
simulations. In the following these galaxy catalogues are 
called input catalogues. 

The input catalogues are used to create synthetic im- 
ages using the programme skymaker^ b y E. B erlin. A 
short description is given in lErben et alJ l|200l|) . Twelve 
30' X 30' images resulting in a 3 deg^ survey are obtained 
for the twelve catalogues. These images are treated in ex- 
actly the same way as real data (like object detection, 

^ Available at: 
http: / / terapix.iap.fr / cplt / oldSite / soft /stuff/ 

® Available at: 
http: / / terapix.iap.fr/ cplt / oldSite/soft /skymaker/ 



4.1. Noise-free case 

To neglect the intrinsic shape of galaxies we perform a 
cosmic shear analysis of the input catalogue for intrinsi- 
cally round objects using the aperture mass statistics. The 
results are given in Fig. 151 The displayed A CDM predic- 
tion is based on the lPeacock fc Doddsl ^^199fi^ model of the 
non-linear power spectrum. The cosmic shear signal is in 
excellent agreement with the theoretical prediction. Hence 
we conclude that the pipeline measures the two-point cor- 
relation function and calculates the aperture mass disper- 
sion exactly, and does not create artificial B-modes. 

For angular scales smaller than 6q < 2 arcmin the 
E- and B-mode signals are significantly smaller than 
expected. This effect is explained and quantified in 
iRilb inger et alJ l(200fTl) . The measured two-point correla- 
tion functions, are set to zero for small angular scales 
(here we choose the same value as for the observations: 
^0 < 6 arcsec) because close galaxy pairs are rejected. This 
lack of shear correlation measurements on small scales re- 
sults in a mixing of E- and B-modes. 

Note that for angular scales larger than two arcminutes 
the marginal difference between the theoretical prediction 
and simulations can easily arise from cosmic variance due 
to the fact that we only use 3 deg^ from the numerical 
simulations. 

4.2. Realistic case 

Using the "old" pipeline the shear measurements of sim- 
ple skymaker simulations yielded, on average, an un- 
derestimation of__the shear by ~15% (STEP 1, see 
Hevman s et al J 12006a') . With our (modified) pipeline we 
are currently able to recover the shear with less than w 3 % 
calibration bias (depending on the PSF model). This has 
been checked with a second set of more realistic shapelet- 
simulations (STEP 2: Massey et al. in preparation). A 
short comparison between the modified and "old" pipeline 
is given in table El 

The same analysis as in Sect. I4.ll is performed for the 
output catalogue. We present the results of the analysis 
done with the "old" and modified pipeline, prior and after 
the changes we have made to our pipeline due to STEP 1 




Fig. 3. Noise-free case: cosmic shear analysis of the input 
catalogue. The number density of background sources is 
n = 15arcmin~^, and the sources are intrinsically round. 
The signal in each bin is the average of the 12 fields, where 
the statistical weight for each field is unity; the error bars 
are obtained by bootstrapping, see Sect. l6.ll and are highly 
correlated. Displayed is the aperture mass dispersion (E- 
and B-modes). The line is a ACDM prediction assuming 
rim = 0.3, = 0.7, r = 0.21 and erg = 0.9 and Zsourco = 1 
(Peacock & Dodds model of the power spectrum). 

Table 2. Comparison between the modified and "old" 
pipeline. Both pipelines are based on the KSB algorithm 
outlined in the text. The c hanges made to the pip eline 
are based on results given in lHevmans et al.l l)2006a(l and 
Massey et al. (in preparation). Key: this change was 
made after STEP 2. 



modified pipeline 


"old" pipeline 


anisotropy measurement^: 




q* matched with a Gaussian 


fixed Gaussian filter 


filter scaled to the size of the 


scale 


galaxy image to be corrected 




cuts: 




l7l < 1-0 


I7I < 0.8 


tj:[P^] > 0.05 


trip*^] > 


calibration factor 7corr ~ 7/O.88 


no calibration factor 



and 2. Both results are compared with the expected signal 
of the noise-free case, see Fig. 0] Clearly visible is the im- 
provement of the modified pipeline compared to the "old" 
pipeline. However, we still slightly underestimate the cos- 
mic shear signal with the modified pipeline by a few per- 
cent. In addition, our pipeline does not create significant 
artificial B-modes (note that the error bars in Fig. 0] are 
strongly correlated and the B-mode signal varies around 
zero). 

We conclude that with our updated cosmic shear 
pipeline including shear estimation and calculation of var- 
ious cosmic shear statistics we can perform accurate mea- 
surements of the cosmic shear signal within the currently 
attainable accuracy. 



Fig. 4. Realistic case: cosmic shear analysis of the output 
catalogue. The same objects are used as in Fig. El Shaded 
area: 1 a range of the signal in Fig. O Displayed are the 
aperture mass dispersion (E- and B-modes) and (M^). 
In addition, the Eoid-mode signal obtained via the "old" 
pipeline prior STEP 1 is displayed. 

5. Dealing with systematics in the GaBoDS data 

In this section we describe various tests to discover system- 
atic errors in our data set. First we present four simple and 
quick quality checks of the data. Second, we present a pow- 
erful standard test for systematics, the cross-correlation 
bet ween uncorrected st ars and anisotropy corrected galax- 
ies (|Bacon et al.ll200,sj) . All tests described basically pro- 
vide checks on the quality of the PSF anisotropy correc- 
tion. However, to test the smearing correction (i.e. the 
correction for the PSF size) we have to rely on simula- 
tions. For that purpose the aforementioned STEP project 
has been initiated. 

5.1. Four basic quality tests 

I. Shear - PSF anisotropy. We calculate the average 

1 ^ 

K2 7i,2) = ^E^tr'(^O7i,2(0.) (29) 

1=1 

before and after anisotropy correction. The quantities 71,2 
are the shear components of the galaxies and e*''2°^ are the 
components of the polynomial fit function (see Sect. I3.4|l 
at the position 6i of the galaxy. The result is displayed 
in Fig. [S| For some data sets the average, {e*-fi), is not 
consistent with zero after anisotropy correction. However, 
the average is reduced by at least a factor of 20 by applying 
the anisotropy correction. 

II. Average shear. We calculate the average shear (71,2) 
of all galaxies in the magnitude bin R € [21.5, 24.5] for 
each GaBoDS-field, before and after anisotropy correction, 
see Fig. |5| For most of the fields the average shear val- 
ues, (71,2}, are significantly different from zero before the 
anisotropy correction and 71 is on average significantly 
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Fig. 5. Test of the anisotropy correction in the magnitude 
range R G [21.5,24.5] using Eq. it^ . Left panel: before 
anisotropy correction; right panel: after anisotropy cor- 
rection. Note the drastic change of scale. 
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Fig. 6. Average shear of all galaxies in the considered 
magnitude bin R G [21.5, 24.5] calculated for all GaBoDS 
fields before (left panel) and after anisotropy correction 
(right panel). Because of clarity we only plot one error 
bar for a single field (thin lines). The average of all fields 
is plotted with a solid error bar. Note the drastic change 
of scale. 



larger than zero. The scatter is strongly reduced after 
anisotropy correction. For most of the fields the average 
shear is consistent with zero. 
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Fig. 7. Averaged shear components versus the ellipticity 
components of the PSF at the galaxy position (polynomial 
value) with and without anisotropy correction. Upper 
two panels: AMl-field; lower two panels: C04pl-field. 
The C04pl-field shows still a small bias (71 and 72 are on 
average larger than zero) after the anisotropy correction. 



III. Average shear in bins. A further simple test of resid- 
ual systematics is given in Fig. [T] where we display for 
two fields the average shear as a function of the PSF- 
anisotropy (value of the polynomial fit function at the 
galaxy position). The two examples in Fig. show that 
the anisotropy correction works well for the AMl-ficld, 
whereas the C04pl-field shows still a small bias in both 
shear components (the average shear values are on aver- 
age larger than zero) after the anisotropy correction (see 
table lA.ll for details of individual fields). 

IV. E- and B-modes of stars. In order to estimate the 
amount of residual systematics in the data set, we utilise 
the stars used for the anisotropy correction. Therefore the 
correlation functions ^± are measured for all stars before 
and after the anisotropy correction and the correspond- 
ing E- and B-modes ((Mfp) and {Mj_)) are calculated. 
The results are presented in Fig. |H| After the anisotropy 



correction we still measure a significant signal which is, 
however, ten times lower than before correction and ten 
times lower than the expected cosmic shear signal. 

Notable is the significantly larger E- than B-mode sig- 
nal of uncorrected stars indicating that on average the 
PSF-anisotropy pattern influences the E-mode signal of 
galaxies more strongly than the B-mode signal. 

5.2. Cross-correlation of galaxies and stars 

5.2.1. ^gyS 

A way to estimate the amount of systematics left in the 
cosmic shear signal (like an insufficient PSF anisotropy 
correction) is given by the cross-correlation of uncor- 
rected ellipt icities of stars and PSF-corrected ellip ticities 
of galaxies l|Brown et al.ll200,'i iBacon et al.ll200,4 . Since 
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Fig. 8. E and B-modes of uncorrected (upper panel) and 
anisotropy corrected stars (lower panel). The signal of 
the corrected stars is significantly inconsistent with zero. 
Note, however, the cosmic shear signal is about 10 times 
larger. 

we use the correlation functions ^± to estimate cosmolog- 
ical parameters we propose a similar estimator: 



?SYS (^) 



Nt 

1=1 



*7) + (g)l(e*7) + (^)l 



where TVf is the total number of fields, (|e*l^)i is the 
auto-correlation function of stars and Wi is the statisti- 
cal weight of a single field obtained by bootstrapping on 
galaxy basis, see Sect. 16. ll The cross-correlation function 
is given by 



E5 «'«,T^.,±4,.7x.,)A.,(e) 



(31) 



where we set w* ~ 1 and calculate Wj^ with Eq. if^ . 
The quantity e* is the uncorrected stellar ellipticity and 7 
is the shear estimate of the galaxies. 

The definition of Csys t&kes into account that the 
cosmic shear signal, ^4., can be negative for individual 
fields. In combination with the normalisation of the auto- 
correlation function of stars, ^^ys directly comparable to 
the signal, <^+. The results are displayed in Fig. ^2 The sig- 
nal ^+ is about 10 times larger than the cross-correlation 
signal CsYS- conclude that our way to correct for the 
PSF anisotropy is sufficient to obtain a clean cosmic shear 
signal using the two-point shear correlation function. 



5.2.2. Mcross;E and Mcross;B 

We calculate the two-point cross-correlation function 
(e*7)^ for every field using Eq. (|31|l . The cross-correlation 
function (e*7)^ is inserted in Eq. (|14|l and Eq. H15|) to ob- 
tain 



and 

A/cross;B = (Mx ^^x ) i^) = ^N, /flX 



,(33) 



respectively. The quantity A^f is the total number of fields 
and Wi is the statistical weight of a single field obtained 
by bootstrapping on galaxy basis, see Sect. 16.11 The re- 
sults are displayed in Fig. 1^ and 1101 for the different data 
sets and for the total data set, respectively. The cross- 
correlation signal, Mcross;E, is significantly non-zero for 
the CO-fields and shows the same characteristic as the 
E-mode signal. Hence our anisotropy correction seems to 
work o nly insufficiently for the se fields. Note that, in con- 
trast to lSemboloni et all l(2n0f]|) . a function similar to £,^ys 
(a function normalised with the auto-correlation of stars, 
(|M*p)) is not calculated for Map. This is because in 
our case we analyse single fields which are shallower and 
smaller in size, hence the Map auto-correlation function of 
stars is very noisy and can have zero-crossings. 



6. Cosmic shear analysis 

6.1. Error estimates and field weiglits 

Bootstrapping on galaxy basis. To obtain the statistical 
(30) weights of (Mfp) and £,± for each field i and angular bin j 
we make 200 bootstrap samples of the galaxy catalogue as 
follows. We randomly draw iVgai times galaxies from the 
galaxy catalogue for each field i with putting back (where 
A'gai is the total number of galaxies in the field i) and place 
them at the same position with the same orientation as 



before. For each bootstrap sample, (Afap)ij and £,±,ij are 
calculated for each angular bin j. The statistical error for 
each field and bin is then the bootstrapping variance afy 
The average measurement signal of the total survey in 
each bin j is calculated as 



and 



(34) 



(35) 



where iVf is the total number of fields. As weight Wij we 
use the reciprocal bootstrapping variance Wij = l/crfj- 
In Fig. lA.ll (App. ^ the relative statistical weights of 
the GaBoDS fields calculated for {M^p) versus the angu- 
lar bin (aperture radius Oq) are displayed. The statistical 
weights of fields with a low number density of galaxies 
are lower than of those with a high number density, as 
expected. This can also be seen in Fig. IA.2I (upper and 
lower left panels) where the square root of the variance, 
aij is plotted for two fields with different number densities 
(FDF-field: 1.9 x 10^ galaxies and C04p3-field: 1.3 x 10"* 
galaxies). Clearly visible are also the large error bars for 
small angular scales since the number of pairs is small for 
them. 
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Fig. 9. Measurements of the aperture mass statistics in i?-band for all six data sets. Left panel : E-mode (squares) 
and B-mode (triangles). Right panel: Afcross;E (squares) and Mcross;B (triangles). For clarity E- and B-mode points 
are slightly offset. 



Bootstrapping on field basis. The measurement error of 
the combined signal for each angular bin j is obtained 
by 2000 bootstrap samples of the field sample as follows. 
We randomly draw Nf times fields out of the field sam- 
ple with putting back. This bootstrap sample is combined 
according to Eqs. and With the 2000 bootstrap 
samples we estimate the PDF of the statistical errors in 
the combined signal including cosmic variance and obtain 
covariances of the errors for the final analyses. With the 
PDF we obtain for each angular bin j the error of the com- 
bined signal by calculating the 68% confidence intervals 
about the mean. For equal weights the calculated 1 a vari- 
ance in each bin is comparable to the usual field-to-field 
variance in each bin as the PDF of the bootstrap sample 
is almost symmetric in our case. 

6.2. Cosmic shear signal 

A powerful way to reveal possible systematic errors is the 
application of the aperture mass statistics as it provides an 
unambiguous splitting of E- and B-modes. The presence 
of non-vanishing B-modes is a good indicator for system- 
atics arising, for instance, from an imperfect anisotropy 
correction. However, Kilbinger et al. (2006) found in their 
work a mixing of E- and B-modes if a cut-off in £+ on 
small angular scales occurs. According to iKilbinger et all 
l)2006|) the deviation of the biased E- and B-mode (the 
E- and B-mode signals are smaller due to the cut-off) 
from the true E- and B-mode signal is only f % for an- 
gular scales up to 4 arcmin for a cut-off in ^± at 6 arcsec. 
In our work we calculate {M^ ) from ^± in the inter- 



val [6", 40']. For the parameter estimation we only con- 
sider the E-mode signal from angular scales larger than 
4 arcmin (see below), so the bias described is negligible. 
ISchneider fc KilbingeJ l)2006|) overcome the general prob- 
lem and introduce the ring-statistics which unambiguously 
splits E- and B-modes. 

At this stage it is worth mentioning that B-modes 
can also arise from the i ntrinsic alignment of g alaxies 
(e.g. iHeavens et all l2000t ICrittenden et alJ l200lif which 
could be removed, in principle, bv using redshift informa- 
tion jKing fc Schneideill2002L l200,4 iHevmans & HeavensI 
bOO,*^ . A source of negative B- modes is the intrinsic sh ape- 
shear correlation predicted bv lHirata i^eliak ( 2004) and 
analyse d with numerical simulations by .Hevmans et alJ 
120063). ]n addition, B-modes can arise from clus- 
tering of galaxies fro m which the shear is measured 
l)Schneider et al.ll2002|) . All three sources of contamina- 
tions of the B-mode signal combined are probably not sig- 
nificant as the statistical errors are still quite large. In 
Sect. l6.5l wc discuss the contamination of the E-mode sig- 
nal by the intrinsic alignment and the intrinsic shape-shear 
correlation of galaxies which can be much larger and sig- 
nificant. 

As pointed out in Sect. 01 we split our data set into six 
different sets depending on the data source. The data is 
heterogeneous with respect to the seeing conditions, ob- 
serving strategy and other characteristics (see Sect. 13. ij) . 
Thus splitting of the data into different sets allows us to 
reveal potential systematics depending on the data source. 
In Fig.inithe average E- and B-mode signal and the average 
signal of the cross-correlation between uncorrected stars 
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and PSF corrected galaxies for E- and B-modes, Mcross;E 
and Afcross:B, of the different data sets are displayed. In 
particular the cross-correlation signal, Mci-oss;E: of the CO- 
fields is significantly non-zero and shows the same char- 
acteristic as the E-mode signal. In addition, the B-mode 
signal within the interval € [2', 7'] is significantly larger 
than zero, hence we exclude those fields from our final 
analysis. 

By excluding the CO-fields we obtain, in the magni- 
tude range R G [21.5,24.5], 7.8 x 10^ galaxies, corre- 
sponding to a number density of n = 16 arcmin"^. Taking 
into account the weights of individual galaxies this re- 
sults in an effective number density of source galaxies of 
rioff — 12.5 arcmin"^. In Fig. ^] the average E-mode, B- 
mode, (Mx)-signal and the average signal of the cross- 
correlation between uncorrected stars and PSF corrected 
galaxies of all fields, except the excluded CO-fields, are 
displayed. The average B-mode signal is consistent with 
zero within the 1 cr-range for 6*0 > 4 arcmin, and the cross- 
correlation between uncorrected stars and corrected galax- 
ies, Mcross;B, IS consistcnt with zero. Hence the B-mode 
signal does not suffer from an imperfect anisotropy cor- 
rection. In addition, the (M^)-signal is consistent with 
zero indicating a clean data set, free from any systematic 
errors. The cross-correlation signal, Mcross;E, however, is 
about 1 cr to 2 cr larger than zero for 9q > h arcmin and 3 a 
for 9o G [3', 5']. This suggests that the PSF-anisotropy cor- 
rection is not perfect and biases the E-mode signal even 
though the B-modes are consistent with zero. However, 
Afcross;E and Mcross;B are not normalised with the auto- 
correlation function of stars, so Mcioss;E is not directly 
comparable with the lensing signal itself. Therefore it is 
hard to judge how large the possible impact of systematics 
on the E-mode really is. 

Taking into account the B-niodes, the (M^) and the 
cross-correlation signals Aicross:E and A/cross:B we conclude 
that the influence of systematics on the calculated E-mode 
signal is negligible on angular scales larger than four ar- 
cminutes within the currently attainable accuracy. For 
scales below four arcminutes the B-modes are slightly pos- 
itive (about 1 a) and the Mcross;E-signal has its maximum. 
For this reason, we exclude in our further (M^p)-analysis, 
scales below four arcminutes. 

In Fig. ^2 the average and ^^ys'^ig^i^ls of 

all fields, except the excluded CO-fields, are plotted. The 
^g'yg-signal is at least a factor of 10 smaller than the 
signal so we conclude that the two-point shear correlation 
function is not significantly influenced by an imperfect 
anisotropy correction. 

Both, the aperture mass dispersion (in the interval 
Oq € [4', 20']) and the two-point shear correlation function 
(in the interval G [0'.8, 33']) are used in the following to 
estimate cosmological parameters. 
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Fig. 10. Various measurements of the aperture mass 
statistics of 52 CaBoDS fields (the CO fields and CAPO 
are excluded, see text) in i?-band (« 13deg^). All galax- 
ies in the magnitude range R G [21.5, 24.5] are used for 
the cosmic shear analysis. The total number of galaxies is 
~ 7.8 X 10^ corresponding to an effective number density 
of ricff = 12.5 arcmin^^. Upper panel: E- and B-mode 
decomposition of the Map-signal. The line is a ACDM pre- 
diction assuming Vl^ = 0.3, JIa = 0.7, F = 0.21 and 
CTg = 0.9; the redshift distribution function is given by Eq. 
(|36|l . For clarity the B-mode points are slightly offset to 
the right. Lower panel: Cross-correlation between un- 
corrected stars- and corrected galaxies for E- and B-modes 
and {Ml). For clarity the {Ml)- and (M*pM|^')B-points 
are slightly offset to the left and right, respectively. 



6.3. Photometric redshifts 

The GaBoDS data set includes 2deg^ of deep UBVRI- 
band observations from the DPS, which yields photomet- 
ric redshift information for about 8% of the objects con- 
sidered for the cosmic shear analysis. The redshift cata- 
logue was o btained using the hyper -z public software de- 
veloped bv iBolzonella et al.l l)2000(l . For the photometric 
redshift estimation we used the observed galaxy templates 
by Coleman et al. (1980). 

The synthetic galaxy templates derived from the li- 
brary of lBruzual A. & Chariot) lll993'> yield a smaller scat- 
ter in the difference between photometric and spectro- 
scopic redshifts for galaxies with magnitudes i? < 23. 
However, the redshift distribution of galaxies in the mag- 
nitude interv al R G [21.5,24 . 5], obt ained with the galaxy 
templates bv lColeman et akl ()l980() . matches the redshift 
distribution obtained from the GOODS-MUSIC sample 
(Grazian et al. 2006) slig htly better than that obtain ed 
with galaxy templates bv lBruzual A. fc CharlotI 1 19931) . 
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Fig. 11. Measured two-point correlation functions ^+ and 
^_ . In the lower panel CsVs (^^® -^1- ^'^^^^ plotted addi- 
tionally (triangles). The line is a ACDM prediction (same 
parameters as in Fig. 110(1 . 



For the Chandra Deep Field South (name in our pa- 
per: AXAF), a database of spectroscopically observed 
objects from the VIMOS-VLT-Deep-Survev (VVDS, 
iLe Fevre et"al]l2004 exists. A comparison between the 
redshifts of 500 VVDS-objects and photometric redshifts 
obtained with hyper-z, in the considered magnitude inter- 
val R e [21.5, 24.5], reveals a systematic underestimation 
of the photometric redshifts of 0.055 (1 -I- (z)). This is cor- 
rected for in an ad- hoc manner by adding 0.055 (1 -|- (z)) 
to each photo-z of a galaxy. 

To estimate the redshift distribution of the background 
sources we associate the redshift catalogue with the lens- 
ing catalogue. In this way the weighting of individual 
galaxies (see Eq. I28II is taken into account in the estima- 
tion of the redshift distribution. We only consider those 
galaxies with photometric redshift estimates in the inter- 
val z G [0,3]. In addition, we only take those photo-z es- 
timates into account with a hyper-z 68 %-confidence in- 
terval smaller than 0.55 (1 4- z). With this cut we exclude 
objects with highly uncertain redshift estimates which are 
also galaxies with a low signal-to-noise ratio in i?-band, 
see Fig. [El 

We obtain 6.2 x 10^ objects with accurate photometric 
redshifts corresponding to 70 % of the DPS-lensing objects 
and about 8 % of all lensing objects. With these objects we 
estimate a redshift distribution of the DPS-lensing galax- 
ies in the given magnitude bin. To acquire a smooth red- 
shift distribution for all galaxies of the GaBoDS galaxy 
lensing catalogue we fit the following dis tribution function 
pat, introduced bv lBrodwin et alJ 1 20061) . to the measured 
redshift distribution, 




Fig. 12. Upper panel: Relative hyper-z 68 %- 

confidence interval, conf/(l-|-z), of lensing galaxies and 
its dependence on signal-to-noise ratio (snr) in i?-band. 
Lower panel: Relative confidence interval and its de- 
pendence on photometric redshift z. The vertical line in- 
dicates the cut described in the text. The photometric 
redshifts are obtained from 7 WFI fields of the DPS (with 
5 broadband filters: UBVRT) in the magnitude interval 
R e [21.5,24.5]. 



where H denotes the Heaviside step function. 



and 

P2{z) 



exp 



Pi[zt); 



the normalisation A is obtained by 



dzpi(z)+ / dzp2{z)=^l. 



Pfit(z) =A [pi{z) H{zt - z) + p2{z) H{z - zt)] , 



(36) 



(37) 



(38) 



(39) 



The fit function was binned in the same way as the mea- 
sured distribution. For zt = 1 we minimise ^nd obtain 
as best-fit parameters zq = 0.27 ± 0.14, a = 2.84 ± 0.96, 
/3 = 1.40 ± 0.37, 7 = 2.34 ± 0.53 and zi = 2.16 ± 0.22 
which results in a mean redshift of z — 0.78. The given 
errors are 1 a statistical uncertainties. 

In Fig.^Jthe observed and fitted redshift distributions 
of the DPS galaxy sub-sample in the magnitude interval 
R e [21.5,24.5] are displayed. Note the dip at z = 1.1 
in the observed redshift distribution. This dip is present 
in the estimated redshift distributions of each DPS field, 
yet not for the GOODS-MUSIC sample which comprises 
infrared bands. We assume that the dip is probably a sys- 
tematic effect due to the lack of infrared information. We 
therefore downweighted this redshift bin for the fit. 

A more proper way to determine the redshift distribu- 
tion from the photometric redshift estimates would be to 
deconvolve the latter with the width of the phot-z uncer- 
tainties. However, the corresponding change in the "true" 
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Fig. 13. Normalised number density of galaxies as a func- 
tion of photometric redshift. Solid histogram: measured 
redshifts, shaded histogram: fit given by Eq. The 
photometric redshifts are obtained from 7 WFI fields of 
the DPS (with 5 broadband filters: UBVRI) for 6.2 x lO'^ 
objects (« 70% of the lensing objects) in the magnitude 
interval R G [21.5,24.5]. To better display the tail of the 
distribution we plot inside the diagram a logarithmic plot 
of N{z). The error bars are the field-to-field variance of 
the 7 WFI fields. 

redshift distribution would be fairly small, yielding a cor- 
rection factor which is much smaller than the statistical 
uncertainties of the shear measurements. We have there- 
fore not attempted such a deconvolution and do not con- 
sider this negligible bias of the redshift distribution in the 
following error estimation of cts fSect. 

6.4. Parameter estimation 

The measured cosmic shear signal is related to 
the dark matter density power spectrum, which in 
turn depends on cosmological parameters and the 
source redshift distribution. To estimate the parameters 
OL = (zq, zi, a, /3, 7, CTgi ^^m, h, from our measured data 
vector d (either the binned (Mfp)- or ^±-values) we use 
the generalised likelihood. The probability distribution 
function (PDF) of the model parameters a. given the p- 
dimensional data vector d is called the posterior PDF 
P{a.\d). It is obtained using Bayes theorem as follows 

P{a.\d) cx Pprior(a)/'(d|a), (40) 

where P{d\a) is the usual likelihood function, Pprior(c«) 
is a prior of ol and m{oL) is the p-dimensional model vec- 
tor. Assuming that the noise residual n — d — m obeys 
Gaussian statistics the log-likelihood function reads 

lnP(d|a) = -^(d - mfC'^id - m) + const. (41) 



An unbiased estimator of the inverse covariance matrix is 
given by 

where p is the size of the data vector (in our case 
the number of angular bins, with p = 15), A^f is the 
number of statistical independent data vectors (in our 
case the number of fields, A'f = 52). Without the factor, 
(Aff — p — 2)/{N{ — 1), the error contours would be under- 
estimated and strongly depend on the number of bins, p. 
For a detailed discussion of how to obtain unbiased esti- 
mates of inverse covariance matrices, we refer to Hartlap 
et al. (in preparation). The quantity 

= E (d'f^ " ^"^)(4'^ - d,) (43) 

-'^boot — 1 

is the estimated noise covariance matrix with the angular 
bins i,j G {1, 2, ...,p}. The quantity d is the average signal 
of all fields in the considered bin, see Eqs. (|34() and (|35() . 
and A^boot — 2000 is the number of bootstrap samples on 
field basis (Sect. EHJl. 

In many other works on cosmic shear the covariance 
matrix is computed analytically (especially for those cos- 
mic shear surveys consisting of a large contiguous field). 
The covariance matrix is decomposed as C = Cn -|- Cg, 
where Cn is the statistical noise and Cg is the cosmic vari- 
ance covariance matrix. The cosmic vari ance covariance 
matrix is then computed according to ISchneider et alJ 
(2002) assuming a fiducial cosmological model and an ef- 
fective survey area. Note that in contrast to this approach 
we estimate our covariance matrix directly from the data 
without any further assumptions. This is possible as our 
fields are widely separated in the sky; hence they are sta- 
tistically independent. 

To calculate the model vector m{a.), we have to spec- 
ify the dark matter density power spectrum. Our mea- 
surements probe scales which are affected by the non- 
linear growth of structure. We calculate the non-linear 
power spectrum using t he fitting formula proposed by 
IPeacock fc DodcQ l|l996() . The fitting formula is based on 
a semi-analytical prescription for the evolution of the dark 
matter power spectrum. Although this formula is rela- 
tively simple it reproduces the main features of the stan- 
dard cosmological model and is accurate enough for our 
purposes, considering the large statistical errors of our 
measurements. To estimate cosmological parameters we 
assumed two simple ACDM models: 

— A: a flat universe: + ^^A — 1 and 

— B: a A-universe with JIa, G [0, 1.5]. 

For both models we assumed a negligible baryon content 
(the shape parameter i s given by P = f^m h and the trans- 
fer function is given in iBardeen et al and a primor- 
dial spectral index of n = 1. The remaining free parame- 
ters in our following likelihood analyses are therefore: 
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(44) 



— the five fit parameters of the redshift distribution 
(zq, zi, a, P, 7) with the strong constraint from our ob- 
served redshift distribution, 

— Hubble constant with a strong constraint of 
h = 0.7± 0.07 (supported by the HST key project), 

— mass power spectrum normalisation erg, 

— matter density with the weak constraint fim > 0, 

— in addition for model B: J^a- 

To apply Eq. I|40|) we have to specify the probability func- 
tion of the priors, -Pprior(<^)- It is given by a multivariate 
Gaussian distribution 

Pprior(a) oc exp |-i(a - a)*Cp^[^,.(a - a) 

The correlation between the five fit parameters of the red- 
shift distribution is taken into account in the prior covari- 
ance matrix, Cprior- 

6.4.1. MCMC 

We do not calculate the likelihood function on a grid. It 
would be eight- or nine-dimensional in our case (erg, ilm, h 
or (Tg, fim, h, f^A, depending on model A or B, and the five 
fit parameter of the redshift distribution) and would there- 
fore be too time consuming. The Monte Carlo Markov 
Chain (MCMC) technique overcomes this computational 
limitation. The likelihood function is no longer calculated 
at fixed positions but at randomly selected positions along 
the Markov Chain. 

The MCMC code in our wor k is based on 
the M etrop olis-Hasti n gs alg orit hm, see iMetropolis et al.1 
l|l953(l and iHastingd jigjcil), or i Tereno et all l|2005h for 
recent implementations to estimate cosmological param- 
eters. It performs a random walk through the parame- 



ter space Si = {a]^ 



(i) 
8/9 



) with i G [l,iVs] {Ns is the 



number of accepted steps). The starting point, 5*0, can 
be somewhere in the parameter space. We choose it to 
be close to the expected likelihood maximum. Then a 
candidate for the next step S'i+i, using a proposal PDF 
Q{Si, Si+i), is drawn. The value Si+i is accepted if the 
following inequalities are fulfilled, either 

where a; G [0, 1] is a random number drawn from a uniform 
distribution. If the candidate is accepted we move to S'i+i 
and repeat the procedure. 

In general the choice of Q{Si, Si+i) is free. To achieve 

convergence in a reasonable time we choose, as in 

iTereno et al.l l|2005f) , a multivariate Gaussian distribution 

Q{S^,S^+l) cx exp - 5,+i]*M-i[S, - S,+i]l. (46) 



For M we choose during a pre-phase (first 1000 steps along 
the chain) the Fisher matrix Fij at the position Sq, the 
given starting point: 



_ ^prior 



-I- 



k,l 



doti 



\kl 



drrii 



(47) 



where 



pprior _ 



<~i o prior 

aaiOai 



(48) 



is the Fisher matrix of the priors. After the pre-phase we 
replace M by the covariance matrix, cov(ai,aj), of the 
parameters estimated from the first 1000 steps. As Q in 
Eq. igni satisfies Q{S^,S^+l) = Q{S,+i,Si), Eq. ^ sim- 
plifies to PiS^+i\d)/P{S^\d) > X. 

Usually several independent MCMCs are started. For 
every chain the first 5*^ values are rejected until the 
MCMCs have forgotten about their initial start position. 
In this way a possible bias due to the starting position 
is avoided. To estimate the size of the so-c alled burn-in- 
phase we use the Gelman & Rubin test ( Gelman fc RubinI 
Il992l) . For each parameter, aa, first the variance of the 
parameters, tr^^^ , inside every chain fc, and then the aver- 
age variance (t^^ = 1/fc J^k '^oki all chains is calculated. 
Next the average of each parameter, aki , inside each chain, 
and then the variance of all averages, tr? , is calculated. 
The average, ctki, and variance, cr^^^ , inside each chain 
is calculated using only S'i, with i e [n, 2n]. For n 00 
both variances, CTq^ and (j^^ , converge to the same value. 
Therefore the ratio 



R 



n- 1 



(49) 



calculated for each parameter ai is a reasonable con- 
trol param eter for the convergence of each chain, see 
ITereno et a l. (2005). As long as R is larger than 1.2 we 
reject all Si with i < n. 

To decorrelate the Si values a thinning is applied. 
O nly every fifth value of the chain is used as proposed 
bv ITereno et al.l l|200fJ ^. 

We marginalise over the Hubble parameter h and the 
fit parameter of the redshift distribution, zq, a, (3, 7 and 
zi, meaning that we project all points of the MCMC onto 
the (ri,„,crg), (ilAjCTg) or (ri,„, r2A)-plane. Compared to 
the prior, the Hubble parameter and the redshift dis- 
tribution parameters cannot be improved with the data 
at hand. To determine the {la, 2 a and 3 a) contours, 
the point number density in the (Sim, erg), (f^A,crg) and 
(f2„i, r2A)-pla'ne is smoothed and the logarithm of the 
number density map is calculated. The maximum log- 
density is subtracted everywhere. In this way — Ax^/2 = 
InL — Ininiax is obtained. 

6.4.2. Results 

Estimated parameters. For a flat universe (ACDM model 
A) both the aperture mass dispersion and the two-point 
shear correlation function is used to constrain cosmologi- 
cal parameters. For (M^^) and ^± we ran four independent 
MCMCs resulting in 15000 points (excluding 1500 points 
of the burn-in phase) and 29400 points (excluding 1000 
points of the burn-in phase), respectively. As a result we 
obtain the joint constraints on f^ni and erg, see Fie. 1141 and 
Fig. El The confidence contours of erg and fim reveal the 
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Table 3. Joint constraints on and erg from our cos- 
mic shear analysis using the two-point shear-correlation 
function and the aperture mass dispersion. We assume 
Q\ + = 1, r = fini n — 1, rijjj G [0; 1] ^^id 
marginalise over the Hubble parameter and the redshift 
distribution as described in the text. Stated below are me- 
dians and 1 a errors. 











Qm{as = 0.8) 


(a4> 




^■"^-0.20 


0.93 ±0.14 
0.80 ±0.10 


0.37 ±0.08 
0.31 ±0.07 





Fig. 14. Model A: joint constraints on and as ob- 
tained from our cosmic shear analysis of the GaBoDS data 
set using the two-point shear-correlation function and the 
Peacock & Dodds model for the non-linear power spec- 
trum. We assume flA + = 1, T = ^mh, n — 1, 
£ [0, 1] and marginalise over the Hubble parame- 
ter and the redshift distribution as described in the text. 
Large panel: black contours: cosmic shear analysis using 
the two-point shear-correlation function, light grey con- 
tours: aperture mass dispersion (shown are the la, 2a and 
3a contours). Upper right panel: solid line: PDF of 
for the marginalisation over ag, dashed line: PDF of Jim 
for as — 0.8. Lower right panel: solid line: PDF of as 
for the marginalisation over fim, dashed line: PDF of as 
for rim = 0.3. 



typical "banana" -like shape reflecting the strong degener- 
acy between these two parameters. Due to the degeneracy 
both parameters are poorly constrained without further 
priors, see table O and solid lines in the small panels of 
Fig. [Hand Fig.lT^ 

Assuming, as a strong prior, either = 0.3 or 
as = 0.8, we obtain stronger constrains on as and fini, re- 
spectively, see tableland dashed lines in the small panels 
of Fig. HI and Fig. US 



Fig. 15. Model A: same as in Fig. E]but now for the 

aperture mass dispersion (black contours). Light grey con- 
tours: comparison to the two-point shear-correlation func- 
tion. 



For the ACDM model B the aperture mass dispersion 
is utilised to constrain cosmological parameters. We ran 
eight independent MCMCs resulting in 39000 points (ex- 
cluding 1000 points of the burn-in phase). The result is 
given in tableland displayed in Fig.^| Again, the confi- 
dence contours of as and reveal the typical "banana" - 
like shape. The contours of model B are almost identical 
to model A(A/ap) for ri,„ > 0.3 and differ strongly for 
f2,„ < 0.3 (the contours do not "close" until erg « 1.9). 

Assuming as a strong prior either Sim — 0.3 or as — 0.8 
we obtain stronger constrains on erg and Q^, respectively 
which is not the case for Ha , see table 0] and grey and 
light grey lines in the small panels of Fig. 1161 Assuming as 
a strong prior SI a — 0.7 the constraints do not improve. 
Hence the dependence on SIa is relatively weak, imposing 
only weak constraints. 

However, a joint analysis of cosmic shear and super- 
novae measurements without assuming a flat universe 
(without priors from WMAP) could provide a tight con- 
straint on S7a. Our measurements exclude a S7a larger 
than unity with 3 a and the supernovae measurement from 
iKnoD etal. (2003.) require a dark energy density (S7a > 0) 
with more than 3 a. 



Estimating parameters with (M^p) or ^± ? The difference 
in the constraints on S^m and as, stemming from the ap- 
plied aperture mass dispersion (Af^p) and two-point shear 
correlation function ^± , can arise from the fact that both 
methods probe different parts of the mass power spectrum. 
Another reason for the difference arises from the fact that 
we fit (Mfp) in the interval 6*0 G [4', 20'] and (± in the 
interval 6*0 e [0'.8,33']. 

The aperture mass dispersion unambiguously sepa- 
rates E- and B-modes in the considered range, 9o G [4', 20'] 
hence this method guarantees the revelation of potential 
systematics in contrast to As our data is B-mode 
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Table 4. Joint constraints on fim, erg and fl^ from our 
cosmic shear analysis using the aperture mass dispersion. 
We assume fl^^ € [0,1.5], T = fli^h, n = 1 and 
marginalise over the Hubble parameter and the redshift 
distribution as described in the text. Stated below are me- 
dians and 1 a errors. 



6.5. Systematic errors 

iBernardea^i et all l(l997f ) showed that for a power law 
power spectrum the cosmological parameters ag and 
are related to the top-hat shear variance (7^) (see Eq. (HI 
and the mean source redshift as: 



add. prior 




0-8 




none 
f^m = 0.3 
as = 0.8 
= 0.7 


29+° '"' 



n 70+0.62 

"•'^-0.36 
70+0.09 
' ="-0.08 

'"-0.36 


U.OO_Q 23 

n QQ+0.33 
u.oo_Q 22 

r\ 07+0.35 
-0.26 




0.5 1 1.5 2 



Fig. 16. Model B: bigger panels on the left: joint con- 
straints on r^in, cTg and Q\ using the aperture mass disper- 
sion for the A-universe: ^lA, G [0,1.5]. Upper right 
panel: black line: PDF of flm when marginalising over 
(J8 and Ha. Grey hne: PDF of fi,„ for JIa = 0.7. Dotted 
line: PDF of fl^ for cg — 0.8. Centre right panel: black 
line: PDF of J^a when marginalising over ag and Om- Grey 
hne: PDF of for f^m = 0.3. Dotted hne: PDF of 
for CTg = 0.8. Lower right panel: black line: PDF of erg 
when marginalising over and ^a- Grey line: PDF of 
CTg for Ha = 0.7. Dotted line: PDF of erg for iln, = 0.3. 
Small panel in the lov^rer middle: comparison to model 
A (flat universe): joint constraints on fim, erg. 



free for Oq > A arcmin, the influence of systematics on the 
calculated E-mode signal is negligible within the interval 
9o G [4', 20']. We conclude that our cosmological param- 
eter estimation using (Mfp) is therefore preferable to 



(7^ 



-1.7 



(50) 



In the following we utilise this formula to obtain a rough 
systematic error estimate of erg . 

Redshift distribution. The redshift of source galaxies has 
a large impact on the estimation of cosmological parame- 
ters, see Eq. (jnOJ. It is therefor e crucial to know the red- 
shift distribution acc u rately . See lHuterer et al.1 l)200(il) and 
I Van Waerbeke et"al] l(200(i^ for recent studies of photo- 
metric redshift errors and their influence on cosmic shear 
analyses. 

Unlike most of the other cosmic shear surveys we es- 
timate our redshift distribution from a sub-sample of the 
total lensing galaxy catalogue and not from an external 
redshift sample. This has two advantages. First, a redshift 
distribution of pure lensing galaxies is obtained. Second, 
the sample variance can be estimated from the field-to- 
field variance of the redshift distributions obtained from 
single fields and can be taken into account in the param- 
eter fit. 

The sources of error of the redshift distribution are: a) 
the sample variance, as we only have redshift information 
available for 8 % of our lensing galaxies, and b) the uncer- 
tain redshift information especially for those galaxies with 
i? > 23 mag. 

To obtain a more reliable redshift distribution, about 
7 % of the faint galaxies with highly uncertain redshifts are 
rejected using the aforementioned filter according to which 
the hyper-z 68 %-confidence interval should be smaller 
than 0.55 (1 -|- z) (note that the hyper-z-uncertainties 
do not follow a Gaussian distribution). With this filter 
the mean redshift is 8 % lower than the mean redshift of 
the unfiltered redshift distribution, since most of the re- 
jected redshifts lie in the interval [0.8,2.4], see Fig. 1121 
However, as there are no spectroscopic or deep photomet- 
ric redshift information obtained from bands which also 
comprise the infrared available it is impossible to judge 
if the rejected galaxies are really high-redshift galaxies 
or, for instance, dwarf galaxies at low redshift. Thus it is 
not possible to say if we underestimate the mean redshift 
of our redshift distribution applying the aforementioned 
filter . Comparing our mea sured mean redshift with that 
from lBrodwin et al.l l)2006l) obtained from galaxies in the 
same magnitude range we are confident that at least the 
mean redshift of our estimated redshift distribution is not 
strongly biased. 

To obtain a rough estimate of the systematic error in 
CTg due to the uncertain redshift distribution we simply 
assume the error of the mean redshift to be of the order of 
8%. This results approximately in an bias of 5% for erg, 
see Eq. (|5n)l . 
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Intrinsic alignment &l intrinsic shape-shear correlation. 

The intrinsic ahgnment of close galaxy pairs introduces 
only a small contamination of the cos mic shear signal 
on sm all angular scales (9 < l.Sarcmin) l)Hevmans et al.l 
The contamination could be removed totally 
by using redshift informa t ion llKing fc Schneider! 120031 
IHevmans fc Heavens! EoO^ . As we do not have multi- 
colour information for most of our fields we cannot correct 
properly for this effect. However, the intrinsic alignment 
contamination is probably negligible as we average over a 
wide range in redshift and only probe angular scales with 
9 > 4 arcmin (in the case of the Map-statistics) . 

A further source of bias is the intrinsic sh ape-shear cor- 
relation predicted b v Hirata fc Seliak! l)2004!) , measured by 
!Mandelbaum et iil] I200(i!^ and analyse d with numerical 
simulat ions bv !Hevmans et al.! l|2006b(l . !Hevmans et al.! 
l)2006b(l predicted that for a survey depth comparable to 
ours, the intrinsic shape-shear correlation produces a small 
negative B-mode and reduces the cosmic shear signal of 
the order of 15% on scales below 20 arcmin resulting in 
an underestimation of erg by about 7%. 

Model. Our estimation of cosmological parame- 
ters is based o n the fi tting formula p ropose d by 
!Peacock fc Dodds! l|l99f]|) . !Hoekstra etal] l|2005!) per- 
forme d a cosmic s h ear an alysis and obtained a trs-value 
using !Smith et alJ l)2003j) which is about 3.5 % smaller 
compared to the erg-value obtained usi ng Peacock fc 
Dodds . As our survey is comparable with !Hoekstra et alJ 
(2005) we assume that we overestimate erg by about 4 % 
using Peacock fc Dodds (assuming the Smith et al. model 
to be closer to reality). 

PSF correction. The presence of non-vanishing B-modcs 
is a g ood indicator for systematics. However, Hockstra 
1)2004!) pointed out that the influence of systematics on 
the cosmic shear signal (E-mode) is not necessarily the 
same as on the B-modes. This is probably the case in our 
analysis as indicated in Fig. |S| The PSF-anisotropy pat- 
tern of the WFI-fields have on average a larger impact on 
stellar E-modes than on stellar B-modes. In addition, we 
saw in Sect. 16.21 that the B-mode signal is zero and does 
not suffer from an imperfect anisotropy correction but the 
E-mode signal probably does. We concluded that the sys- 
tematic error of the E-mode signal for angular scales larger 
than four arcminutes due to an imperfect anisotropy cor- 
rection is in any case significantly lower compared to the 
statistical errors. 

The results in STEP 2 indicate that on average we 
overestimate the shear by 3 % resulting in an overestima- 
tion of the E-mode signal. In contrast, our cosmic shear 
analysis of synthetic skymaker-images in Sect. 0] suggests 
that we underestimate the cosmic shear signal by a few 
percent. Both, the probably small bias of the anisotropy 
correction and the bias in the shear estimation are simply 
accounted for by increasing the systematic error of cg by 
±5%. 



That the average B-mode signal is consistent with zero 
on scales larger than four arcminutes does not imply that 
single fields are B-mode free. It is possible that B-modes 
resulting from different systematic errors only average out 
to zero. If this is the case, the average E-mode value would 
be the same as without any systematic errors assuming 
that systematic errors bias the E-modes in the same way 
as B-modes. However, the full potential of the cosmic shear 
analysis is not achieved in this case as the errors would 
be larger than without any systematics. In our measure- 
ments, the error bars of the average B-mode signal for 
scales larger than 6 arcminutes are smaller than those of 
the E-mode signal. This is expected since B-modes are 
not affected by cosmic variance. Hence, single fields are 
not dominated by systematic effects and we do not in- 
clude this possible increase of the systematic errors into 
the final error estimation. 



Field selection. Our data set comprises one field with a 
massive cluster (A901) which was known before selecting 
the fields. In addition, a field like A901 is unlikely to be 
part of a randomly chosen field sample compared to ours. 
This can be construed as a prior knowledge which intro- 
duces a bias on the cosmic shear signal. Replacing the 
E-mode signal of this field by the average E-mode signal, 
the final signal decreases by roughly 3.5 %. This results in 
an potential overestimation of erg by approximately 2%. 

Error summary. Taking all sources of systematic errors 
into account we see that we already entered a new phase 
of cosmic shear surveys. That is, the systematic errors 
are of the same order of magnitude as statistical errors. 
We obtain for the mass power spectrum normalisation 
assuming fim = 0.3 and including the statistical error 
and the rather conservative estimates of systematic er- 
rors: erg = 0.80(±0.10 ± 0.04 ± 0.04 + 0.06 - 0.03 - 0.01) 
(Icr statistical error, redshift bias, PSF correction bias, 
intrinsic shape-shear correlation bias, model bias, field se- 
lection bias). All given estimates of the systematic errors 
are maximum systematic errors. 



7. Summary and outlook 

We have performed a cosmic shear analysis based on 
15 deg^ of deep high quality i?-band imaging data from the 
Garching-Bonn Deep Survey. Our cosmic shear measure- 
ment is B-mode free within the statistical errors, hence 
there are no significant systematic errors resulting from 
the data treatment (PSF correction, galaxy selection) left 
in the data. This encourages us to use the GaBoDS data 
for further cosmological analyses, such as the galaxy bias 
studies carried out by Simon et al. (in preparation). 

The measured redshift distribution is obtained from 
lensing galaxies from 2deg^ of UBVRI-hsind data of the 
Deep Public Survey, a sub-sample of GaBoDS and not 
from an external redshift sample. 
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The measured redshift distribution in combination 
with the cosmic shear signal is used to perform an un- 
biased estimate of cosmological parameters where, in con- 
trast to many other cosmic shear analyses, we have esti- 
mated the covariance matrix directly from the data with- 
out any further assumptions and have not computed it 
analytically. Assuming a flat ACDM universe with neg- 
ligible baryon content we have derived the mass power 
spectrum normalisation and the matter density while 
marginalising over the uncertainties in the Hubble param- 
eter and the source redshift distribution. As a result we 
obtain erg = OMtoil and n,^ = OAetgi" ^om the Map- 
statistics using the Peacock & Dodds model of the non- 
linear mass power spectrum. For a fixed matter density 
of r^ni = 0.3, we obtain erg = 0.80 ± 0.10 (1 a statistical 
error). Within the error bars this is consistent with re- 
cent cosmic shear measurements, galaxy clusters and the 
WMAP three year result, see Fig. 1171 and table [S] 

We discussed various systematic errors and roughly 
estimated their magnitude. With respect to the magni- 
tude of systematic errors, the most uncertain sources are 
the redshift distribution and the intrinsic shape-shear cor- 
relation. Although the intrinsic shape-shear correlation 
has been analysed with Ai"-body simulations using simple 
galaxy models, the impact on deep cosmic shear measure- 
ments is still unclear. With accurate photometric redshifts 
at hand, measurements of cross-correlation tomography 
could be carried out to provide a useful diagnostic tool 
llHirata fc SelialJ2004l) . To reach the full potential of weak 
gravitational lensing measurements and high accuracy in 
the determination of cosmological parameters, it is there- 
fore essential to have precise redshift estimates to obtain: 

a) unbiased redshift distributions of lensing galaxies and 

b) estimates of the intrinsic shape-shear correlation. 

Our cosmic shear result and the ongoing improvements 
in redshift and shear estimates (e.g. STEP) are quite en- 
couraging for the upcoming wide-field multi-colour sur- 
veys as the future Kilo Degree Survey (KIDS), starting 
beginning 2007. The KIDS aims to image a contiguous 
area of 1500 deg^ in five colours. In combination with 
UKIDS the data will yield photometric redshifts accurate 
to Az/ (1+z) = 0.03 for typical r' — 23.5 galaxies, rising to 
10 % uncertainty at r' = 25. In the primary lensing colour 
the limiting magnitude is predicted to be 7'vcga — 24.3 
(10 (T sky level measured in a circular aperture of 2" ra- 
dius) and the median seeing is predicted to be 0'.'6. The 
depth of KIDS will yield about 10^ galaxies to a median 
redshift of z w 0.8. 

With the KIDS data set at hand we will be able to mea- 
sure the angular power spectrum on large scales. The sta- 
tistical errors of cosmological parameter estimations will 
be reduced at least by a factor of ten, not only because 
the number of background galaxies is larger but also the 
number of galaxy pairs will increase dramatically, in par- 
ticular for large angular scales. With such a large data 
set one can start to probe the equation of state of dark 
energy. 



1.2 rr 




2002 2004 2006 

year 



Fig. 17. Recent determinations of CTg(riin = 0.3) using 
galaxy clusters (triangles) and cosmic shear (squares) in 
comparison to our results (Afapi solid hexagon with bold 
error bars, hexagon with light error bars) and the 
WMAP three year result (pentagon). The open triangle 
and open star on the right are the average of all erg deter- 
minations between 2002 and 2006 (indicated by the ver- 
tical lines) using clusters and cosmic shear, respectively 
(error bars of single measurements are not taken into ac- 
count). The WMAP result of erg would be larger if fim is 
fixed to 0.3. Within a year the measurement points are not 
in chronological order. Values taken from table |S1 Points 
from right to the left in the diagram are associated with 
values from top down in the table. 
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Table 5. Recent determinations of erg using clusters, cos- 
mic shear and WMAP. The WMAP result o f ur would be 
larger if is fixed to 0.3. The result from IChang et al"! 
l)2004) is obtained from a radio survey assuming that the 
median redshift is Zm — 2.2. 



Reference 


0"8 


Method 


our measurements 


0.80 ± 0.10 


WL 






0.93 ± 0.14 


WL 




hembolom et al. (zDOnl 


0.85 ± O.Oo 


WL 


e 


Hoekstra et al. ( 20051 


0.85 ± 0.05 


WL 


e 


Massev et al. (iUUoi 


1 no 1 n 1 c 


WL 


?± 


Jarvis et al. f2006) 


^ — o+O 08 


WL 




Hevmans et al. (20051 


0.68 ± 0.13 


WL 




\T TXT "Ui J- 1 /oriricr\ 

Van Waerbeke et al. (2005) 


0.83 ± 0.07 


WL 


e 


Hevmans et al. (2004) 


0.67 ± 0.10 


WL 




Rnodes et al. (2()()4] 


1.U2 ± 0.16 


WL 




Cnane et al. (2004) 


1.0 ± 0.2 


WL 




Bacon et al. (2003) 


0.97 ± 0.13 


WL 




Brown et al. (2003) 


0.72 ± 0.09 


WL 




Hamana et al. ('20031 


0.78+^.5^ 


WL 


(Mfp) 


Jarvis et al. (zUUo) 


n -71 +0.06 


WL 




Hoekstra et al. (2002) 


n oc+0.09 


WL 


(Mfp) 


Hoekstra et al. (2002b) 


n oi +0.07 

0-81Io.io 


WL 




Van Waerbeke et al. (2002) 


0.98 ± 0.06 


WL 


(Mfp) 


Ketreeier et al. (2002) 


(J. 94 ± (J.z4 


WL 


Rhodes et al, (2001) 


^ „ T _i_n 25 


WL 




^Moli ot al. (2001) 


1.04 ± 0.05 


WL 




Henry (20041 


0.62 ± 0.04 


X-ray clusters 


Allen et al. (20031 


0.70 ± 0.04 


X-ray clusters 


Bahcall et al. (20031 


0.72 ± 0.06 


X-ray clusters 


bcnuecker et al. (2003) 


c\ ''7 a 1 n m 

0.76 ± 0.01 


X-ray clusters 


Piprnnnii pt al fPfin'^l 


n 77 + n (15 


X-ray clusters 


PierDaoli et al. (20031 


0.78 ± 0.06 


X-ray clusters 


Viana et al. (20031 


U-84_o.o3 


X-ray clusters 


Viana et al. (20021 


0.61 ±0.05 


X-ray clusters 


Reiorich & Bohrineer (20021 


0.68 ± 0.04 


X-ray clusters 


Rosati et al. (20021 


0.72 ± 0.02 


X-ray clusters 


Seliak (2002J 


0.76 ± 0.06 


X-ray clusters 


Wu (2001) 


0.91 ±0.11 


X-ray clusters 


Pieroaoli et al. (20011 


1.02 ±0.07 


X-ray clusters 


gukbir & Arnaud (2001) 


0.91 


X-ray clusters 


Blanchard ct al. (2000) 


0.96 


X-ray clusters 


SDergeletal. (20061 


n 74.+o.ob 
'J- '^-0.06 


WMAP 
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Fig. A.l. Relative statistical weights of the GaBoDS fields 
versus the aperture radius obtained via bootstrapping 
on galaxy basis. The COMBO field AXAF is chosen as a 
reference field. The type of each line is the same as in Fig. 

Appendix A: Single fields 

In this appendix we give the main characteristics of the 
single fields and display the relevant Map-statistics and 
relative weighting factors of all fields. 
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Table A.l. List of all fields of the GaBoDS survey. The quantity N is the number of galaxies per field of the final 
catalogue in the magnitude range R € [21.5, 24.5] and 2 rh is twice the measured half light radius of stars in arcseconds. 



name 


KA (J 2000) 


TMr/^ ( toaaa\ 

DhjC (JzUUUj 


tobs[KSecJ 


JN 


2rh[ J 


data source 


code 


A901 


149.07771 


— 10.02734 


18.1 


19513 


0.75 


COMBO-17 


COMBO 


AXAF 


53.13344 


—27.82255 


57.2 


19639 


0.90 


COMBO-17 


COMBO 


FDF 


16.44541 


-25.85742 


11.8 


19169 


0.86 


COMBO-17 


COMBO 


Sll 


175.74860 


— 1.73458 


21.5 


18855 


0.84 


COMBO-17 


COMBO 


SGP 


11.49852 


—29.61047 


20.0 


21465 


0.85 


COMBO-17 


COMBO 


CAPO_DF 


186.03787 


— 13.10764 


13.0 


16014 


0.97 


ESO archive 


OWN 


NDF 


181.36237 


-7.65226 


21.8 


14146 


0.93 


ME IR group 


OWN 


SHARC2 


76.3333 


-28.81805 


11.4 


19101 


0.86 


own obs. 


OWN 


F17_P1 


216.41916 


-34.69460 


10.0 


10345 


1.13 


own obs. 


OWN 


F17_P3 


217.02611 


—34.69463 


10.0 


14135 


0.77 


own obs. 


OWN 




175.25702 


O C ^1 A'y A 

—25.51474 


13.5 


10330 


0.76 


own obs. 


OWN 


A1347_P2 


175.79254 


—25.50918 


7.5 


13738 


0.91 


own obs. 


OWN 


A1347_P3 


175.23976 


—25.00933 


7.0 


15581 


0.89 


own obs. 


OWN 


A1347_P4 


175.79459 


—24.99836 


8.0 


13665 


1.07 


own obs. 


OWN 


F4_P1 


321.65611 


—40.25193 


9.5 


16135 


0.93 


own obs. 


OWN 


F4_P2 


321.71942 


-39.76761 


7.0 


14103 


1.06 


own obs. 


OWN 


F4_P3 


322.32389 


—39.72689 


10.0 


16348 


0.87 


own obs. 


OWN 


F4_P4 


322.32389 


—39.72689 


7.5 


12943 


1.10 


own obs. 


OWN 


Deepla 


343.79506 


-40.19886 


7.2 


16131 


1.00 


EIS 


DEEP 


Deeplb 


343.058572 


-40.22481 


3.9 


9806 


1.30 


EIS 


DEEP 


Deeplc 


342.328125 


-40.20702 


3.9 


13354 


1.03 


EIS 


DEEP 


Deeple 


341.96679 


-39.52874 


9.0 


13374 


1.05 


EIS 


DEEP 


Deep2a 


54.372223 


-27.81551 


6.0 


14785 


1.00 


EIS 


DEEP 


Deep2b 


53.746626 


0*7 0A0/?0 

—27.80862 


5.1 


10177 


1.20 


Elb 


UEEF 


Deep2d 


52.506344 


—27.81774 


3.0 


9745 


1.20 


EIS 


DEEP 


DEEP2e 


53.12291 


—27.30467 


7.5 


16768 


0.97 


own obs. 


DEEP 


DEEP2i 


53.66995 


—27.32400 


7.0 


14494 


1.10 


own obs. 


T^TTITTIT^ 

DEEP 


Deep3a 


171.24559 


—21.68289 


7.2 


14130 


0.91 


EIS 


DEEP 


Deep3b 


170.66159 


—21.70969 


9.3 


14905 


0.92 


EIS 


T^TTITTTT^ 

DEEP 


Deep3c 


170.01909 


—21.69960 


9.0 


15606 


0.90 


EIS 


DEEP 


Deep3d 


169.42875 


-25.85742 


9.3 


14134 


0.80 


own obs. 


DEEP 


AMI 


58.81181 


-49.667762 


7.5 


14956 


0.99 


ASTROVIRTEL 


B8 


B8p0 


340.34886 


-9.59009 


7.2 


13299 


0.87 


ASTROVIRTEL 


B8 


-DOpi 


340.346051 


—9.08957 


4.5 


13654 


0.87 


1 KU V iKl liL 


B8 


B8p2 


340.345309 


—8.58963 


5.4 


13308 


0.87 


ASTROVIRTEL 


B8 


B8p3 


340.345149 


—8.08951 


5.4 
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0.89 


ASTROVIRTEL 


B8 


BSml 
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4.5 
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0.90 


ASTROVIRTEL 


B8 
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<^ A r\ ^ A a A a 
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Fig. A. 2. In each diagram various Map statistics are displayed for each field: cross-correlation between uncorrected 
stars and anisotropy corrected galaxies, {M*M^), cross-correlation between the tangential component and the cross 
component of the shear, Mx , E and B-mode. The name codes are given in table lA.ll For two fields we plotted the 
statistical errors: FDF (upper left panel) and C04p3 (lower left panel). The type of each line is the same as in Fig. 

EH 



